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13 — abstract 

A  comparative  survey  is  made  of  several  estimators  of  a  scale  parameter  of 
the  Beta-Stacy  distribution.  In  essence,  some  optimum  technique  is  sought  to 
permit  more  insight  into  the  statistics  cf  a  joint  distribution  of  two  random 
variables  when  one  set  has  a  generalized  gamma  distribution.  Comparisons  are 
made  on  the  accuracy  lost  when  no  additional  statistics  are  available  on  the 
set  having  the  generalized  gamma  distribution.  This  work  is  a  logical  extension 
of  the  bivariate  warning-time/failure  time  distribution  effort  conducted  by 
Mihran  and  Hultquist  ( i96 7 ) ^ 


Comparisons  among  Estimators  of  a  Scale 
Parameter  of  the  Beta-Stacy  Distribution 
By  N-L.  Johnson*,  TJniversity  of  North  Carolina  at  Chapel  Hill 
and  S.  Kbtz**,  Temple  University,  Philadelphia 


The  Beta-Stacy  distribution,  as  defined  by  Mihram  and  Hultquist 
Cl]  is  the  joint  distribution  of  two  randan  variables,  and  X^,  when 
(i)  X  '  has  a  Stacy  (generalized  gamma)  distribution  with  density 

jL  • 

function 


CD-  (X).  *  [r(a)f1ca“fxiao“1  exp  [-(X]/a)C]  * 


(x^  >  0;  a,c,a  >  0) 


and 


(ii)  the  conditional  distribution  of  given  X^,  is  beta,  with 
parameters  ei>e2  ^  iari ge  0  to  X^,  so  that 


(2) 


6“1  q  ^  ' 

fx  ix  (*iiv =  u-vv  r  -r1 


(°  <  x2  <  x1J  61,  02  >  0) 
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Given  n  independent  pairs  of  random  variables  {  X^.  , 

(j  =  1,2,. . . ,n) ,  each  having  the  Beta-Stacy  distribution  is  define 

It  is  clear  that,  given  {  X. the  observations  {X  }  cm  the  st  ■: end 

2j 


variable  {  X„.  }  provide  no  further  information  on  the  values  of  ’he 
21 

parameters  a,c  and  o  of  the  common  distribution  of  the  X^_.  ’s 


n  c 

In  particular,  if  c  and  a  are  known  then  E  X  .  is  a  suffi  cient 

j=l  ^ 

statistics  for  a,  and  the  minimum  variance  unbiased  estimator  of  a  is 


(3) 


A 

a 


T(na) 


n  c 
(  E  X  C) 


-1  j=l 

T(na  +  c  ) 


Its  variance  is 


(4)  varCa^  =  a2 


I (na)r(na  +  2c"1) 
f:(na  +  c"1)]2 


) 


while  the  Cramer-Rao  lower  bound  for  unbiased  estimators  of  a  is 


(5) 


a2(nac2)-1 


The  ratio  of  var(a^)  to  a2(nac2)~^  tends  to  1  as  n  tends  to  infinity. 
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MLhrara  and  Hultquist  [1]  studied  the  problem  of  estimating  Ct, 

when  observations  are  available  on  {X  .  k  but  nor  on  {X  .  k  it 

23  13 

being  supposed  that  0^  and  0^  are  known,  as  well  as  c  and  a  . 

They  ajiggested  two  estimators,  based  a  the  geometric  and  aiithire'tic 

mean  of  the  X^.’s,  res] actively: 

*3 

rCe^  r  (a)  r(  q1  +  e2  +  n"1) 


(6.1)  sL  = 


-I  n 


r(e  +  n”1)r(a+  (nc)~1r(01  +  02)  J  J 


TT  x  1/n 

«i  23 


.  (0.  +  0  )  r(a) 

C6.2)  a  = - A - 2 - 

5  e1  rCcc  +  c"1) 


1 

n  j 


n 

\  X23 


S,  c  r<0  +  8)r(e,  +  0_) 

Since  E[X  X,]  =  X.5  - i - 1 - SL 

2  1  1 

r(0  +  0  +  s)  r(0  ) 

1  2  1 


and  ECX.^  ]  =  at  r(a  +  sc~^)/T(a)  ^ 


it  follcws  that 

(7)  .  a5+t  T±  *  <Sn)g-1) 


r  (a) 


r(e  +s)r(0n+0  ) 

■ - -  ,  ■  _  r*-  .  *  — , 

r(0i+e2+s)r(01> 


from  (7)  we  obtain 


(8.1)  varCa^)  =  a2  ) 


r(a)  r(0  )r(o+2(ncj  )r(0,r2n  )  ^  rCe  +0 +2n“  ) 

_ 1 _ x _ (  12 

r(e +0.)r(01+Q +2n_1)  |r(a+(nc)"1)r(e  +n~1)  } 2 

l  j 
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(8.2)  «(«,>  =  a^jW!£ 

'(VW 


rCa+ac^irCa) 
CrCo+c-1)]  2 


i 

■) 


A  /A  A  N  2  I  1  +  ^0,)  1 
(8.3)  Cov(a2,a3)  =  a  / - ± _ 


r(a)  r(a+(n+l)(nc)  "S 


1  +  [n>  (§+  fyj]  4  rU+c'1)  r(«+  (n>c)~4^ 


(Formulae  (8.1)  and  (8.3)  disagree  with  the  corresponding  formulae  in 

[1]). 


(9.1) 


(9.3) 


As  n  tends  to  infinity,  we  have 

-2  -2  i 

a  lim  n  var(a  )  =  c  ^(o)  +  ^'(0j  -  ij>'(0  +02) 


-2 

a  lim 
n  -*■  « 


cov(a2,a3)  =  c 


~  "  ^[^(a+c  "S  -  ^'(ot)]  +  020-^  ^0j_+02^  ^ 


where  ^(y)  =  (log  r(y?' 

And,  of  course,  for  all  n 


and  (y)  =  -A-  (i|>(y)). 
dy 


_2  (0-,+l)(b.,  +0.)  r  (a+2c_1)r(a) 

(9.2)  a  n  var  (a  )  =  — - ± - - -  .  - 1 — 2 -  -  1. 

6^0^  0  +1)  {r(a+c  )} 


Table  1  gives  val  *es  of: 


/■»  « 

a""4  n  var(a2),  a“2  n  varCa^  ,  and  corKa^ag)  for  selected  values  of 

c,a,e1,e2  and  n.  The  values  for  n  =  «  are  calculated  from  the  right  hand 

sides  of  (9.1)  -  (9.3).  For  each  set  of  values  of  csa,01  and  e2  ,  n 
is  taken  equal  to  100.  In  a  few  cases  values  are  given  also  for 

A  a  10,  20,  50  and  00  .  These  should  suffice  to  indicate  the 
variation  with  n,  which  is  not  marked. 


In  order  to  see  hew  much  accuracy  has  been  lost  by  ignorance 
of  {X^.  }  the  values  of  a  n  ver  (a.)  (j  =  1,2)  can  be  compared  with 


the  corresponding  value 


-2  .  ir(not)r  (nct!-2c  ) 

a  n  van  (a^  =  n/ 


[r(na+c  ^)]2 


-  1 


for  the  minimum  variance  unbiased  estimator  a^.  These  values  are  given 
in  Table  2. 


Using  the  approximation  \|>(ct)  ==  log(a-  1  )  we  see,  from 

2 

(9*.'l),  that 


(10)  a~2  lim  n  var(a  )  =  c”2(a-  )  -  i|>'(0,+0  ) 

«  2  2  1  ^ 


while 


(11) 


n*  « 


6  — 


£  The  excess  of  varia^)  over  varta^)  ■  can  be  split  into  parts: 

1  ’ 

o"2  {$(<*)  -  a"^}  =  _i_  c^cT^c*-  i  ) 

2  2 

and  ♦,<e.)  -  ♦,(0+Oo) 

The  excess  decreases  as  a  increases  and  as  c  increases;  the  excess  cue  to 
01  and  0j  is  relatively  less  important  when  a  and  e  are  small. 

These  features  can  be  seen  from  the  figures  in  Table  1. 

t- 

Comparison  of  Tables  1  and  2  shows  that  the  variables  {  >  often 

I 

provide  unbiased  estimators  of  c  which  are  of  comparable  acs^ur^cy  Ce.g. 
with  variances  no  more  than  twice  as  great)  to  a-,  ,  Txrr=  is  especially 
notable  for  the  larger  value  of  6n  and  6^  -  another  feature  which  is 
indicated  by  (10)  and  (11).  Except  in  c  cases  (out  of  GO)  in  Table  l, 

1  the  estimator  a2  (based  on  the  arithmetic  mean)  has  a  smaller  variance 
than  ag  (based  on  the  geometric  mean) .  For  the  smaller  values  of 
c  and  a,  the  correlation  between  a„  and  a  is  small  enough  to  give 

^  v 

some  hope  that  the  unbiased  estimator1 


with  A  chosen  to  minimize  var(a  ) ,  suggested  in  [I],  will  gf.  a  use- 

4 

ful  reduction  in  variance.  The  last  two  columns  of  Table  1  »ve  values 
of  A  and  var(a^).  Th.e  reduction  in  variance  is  certainly  worthwhile  for 
c  =  0.5,  but  for  c  =  1.0  and  2.0  is  does  not  seem  to  be  of  much  importance. 
Note  that  as  c  increases  the  value  of  A  neccr.BS  negative. 


-  7  - 


If  6,  and  0  are  large,  with  0,  large  compared  with  0  , 
1  2  .  1  2 

then  X,.  and  X„ .  are  highly  correlated,  and  6,  ^(0.,+0  )X„ .  is  a 
XD  23  1  1  2  2d 

good  unbiased  estimator  of  X^.  .  It  would  seem  likely,  therefore  (in 

view  of  (3))  -fiiat 


-1^-1 


a,  =  (0>0  )  r(no)  [0,  T(nct+  c  )]  (  £  X  .  ) 

^  1  2  1  i=l  23 


would  be  a  good  (though  not  unbiased)  estimator  of  a,  in  such  cases. 
We  note  that  the  statistic 


r(e-1+62+c)  r(e-,) 

r  (e  +c)  r(e  +6  ) 

1  12 


n  o  j-1 


n  c 

£  x„: 

2d 


is  an  unbiased  estimator  of  ac  with  (coefficient  of  variation)  equal  to 


C1W  n-1  /  *  +  1  .  r(y2c)r  (el\  .  x } . 

{  O  CrCSj.+o)]2  r(e1+e2+2c)r(61+02)  ) 


Far  the  minimum  variance  unbiased  estimator  of  a^, 


-1  1  c 

(na)  1  £  X,? 

i=l 


the  (coefficient  of  variation)^  is  (na) .  On  comparison  with  (14),  we 

Q 

see  that  when  c  =  1,  the  efficiency  of  (13) ,  as  an  estimator  of  a  ,  is 

0(1  +  a) 

[1  +  - 2 - ]  . 

h  (91  *  ®2  +  X> 


Table  1  is  based  an  tables  to  5  significant  figures  for 
c  =  0.5,  1.0,  2.0;  a  =  0.5,  1.0,  2.0,  3.0 

»  ®2  =  2.0,  n  =  20>  50,  100,  “ 

calculated  with  an  APL  program  devised  by  Mr.  J.  0.  Kitchen,  to 
whom  we  express  our  gratitude.  Thanks  are  also  due  to  Mrs.  G.  Ballard 
for  assistance  with  the  preparation  of  this  paper. 
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